Quarkyonic Chiral Spirals in the Nambu-Jona-Lasinio Approach 
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We explore the inhomogeneous QCD phases at finite density and temperature using a (3+1)- 
dimensional Nambu-Jona-Lasinio (NJL) model in the large N c limit with an additional attractive 
tensor-tensor interaction channel. For single modulated solutions, the problem reduces to solving the 
gap equation of a Chiral Gross- Neveu (NJL2) theory, whose minimum solution is a periodic array 
of chirally twisted kinks. At zero temperature, the minimum solution reduces to the quarkyonic 
chiral spiral found in quarkyonic matter. The connection between the quarkyonic matter and our 
(3+l)-dimensional NJL model is rooted on the fact that the first reduces to (l+l)-dimensional 
QCD and the second to the NJL2 theory, both of which have continuous chiral symmetry and 
asymptotic freedom. Our findings can be useful to obtain a reliable qualitative picture of the QCD 
phase diagram and the location of the QCD critical point in the technically challenging region of 
intermediate densities and temperatures. 
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Introduction. In this letter we provide a solution to a 
standing problem in the literature related to the existence 
of two different inhomogeneous solutions, each proposed 
as the most energetically favored ground state of large N c 
QCD in the region of intermediate densities and tempera- 
tures. One of these solutions, known as quarkyonic chiral 
spiral [l|, was found in quarkyonic matter |2J]. The sec- 
ond solution was obtained using a (3+l)-dimensional (D) 
NJL theory with conventional scalar and pseudo-scalar 
interactions |3J . The first solution is a complex order pa- 
rameter that forms a chiral spiral of scalar and spin-one 
condensates. The second is a scalar condensate given by 
a real crystalline soliton solution. Both ground states 
are single-modulated and break translational symmetry 
in the direction of the modulation. In the present paper, 
by considering a (3+l)-D NJL model with an additional 
attractive channel that preserves the chiral symmetry of 
the theory, we demonstrate that the most favored ground 
state solution is a crystalline twisted chiral kink that re- 
duces to the quarkyonic chiral spiral solution in the low 
temperature region. This modified NJL model contains 
the NJL theory of [3j as a particular case and is consistent 
with the quarkyonic matter results, suggesting that the 
proposed NJL model can be more physically appropriate 
to describe the intermediate density region of QCD. 

Mapping all the phases of QCD is a goal intensely 
sought after by many theoretical and experimental ef- 
forts [4| . Thanks to the asymptotic freedom of the the- 
ory the most extremes regions of the QCD phase diagram 
in the temperature-density plane are perturbative and 
hence well understood; they are the quark-gluon plasma 
in the high-temperature/low-density corner of the phase 
map and the Color-Flavor-Locked (CFL) superconduct- 
ing phase on the opposite one. At low temperatures and 
densities quarks are confined inside hadrons, whose in- 
teractions can be phenomenologically described by con- 
ventional nuclear physics. 



However, somewhere in the region of intermediate tem- 
peratures and densities, one expects a phase transition 
to occur from a confined to a deconfined phase, where 
gluons and quarks liberate from hadrons. On the other 
hand, confined quarks have a large dynamical mass, while 
in the regions with small coupling the quark mass be- 
comes small. Clearly, by increasing the temperature or 
the chemical potential, a chiral phase transition should 
occur from the phase of heavy quarks to that of light 
quarks. Nevertheless, it is not obvious a priori that the 
confinement and chiral phase transitions have to occur 
at the same temperature [5[. 

Understanding the fundamental physics involved and 
the phases that realize in the strong-coupling region of 
the QCD phase map is nontrivial, as it requires using 
nonperturbative methods and effective theories. To be- 
gin with, one cannot use lattice QCD because of the sign 
problem at finite chemical potential, then the investiga- 
tion of this region can only be done with the help of NJL- 
like [6[ or QCD-like effective models [7J. To pin-down the 
correct physical description of this region, the theoreti- 
cal results will have to be eventually contrasted with the 
outcomes of future experiments at various international 
facilities that plan to explore the intermediate-baryon- 
density and small-temperature region in the near future. 

The strong-coupling region is even more technically 
challenging due to the possibility of inhomogeneous 
phases (IPs). Phases with inhomogeneous condensates 
(IC) have been found in the large N limit of QCD [7[, 
in quarkyonic matter |l| , and in NJL models |3j, I8l4l0 1 . 
These phases seem unavoidable at intermediate densities 
and low temperatures. Indeed, while most NJL mod- 
els predict a first-order chiral transition with increasing 
density |ll[ , it is likely that the transition will occur via 
some intermediate state(s) characterized by inhomoge- 
neous chiral condensates, because the density promotes 
an imbalance between quarks and antiquarks and an en- 



ergy separation between them. Since it costs energy hop- 
ping the antiparticlcs from the Dirac sea to pair with the 
particles at the Fermi sea, the system prefers instead to 
pair quarks and holes giving rise to ICs. 

In this paper we are interested in exploring 1-D modu- 
lated ICs at finite temperature and density in a modified 
NJL model in the large Nc limit that takes into account 
new attractive channels. Previous studies of 1-D mod- 
ulated ICs within (3+l)-D NJL models considered only 
the possibility of chiral and pion condensates [3j, I8l4ld| . 
However, there are other attractive channels that respect 
the original chiral symmetry of the theory and can con- 
dense when some spatial symmetries are either sponta- 
neously or explicitly broken. This is known to occur, for 
instance, in color superconductivity [121 ] and quarkyonic 
matter [13J in the presence of a magnetic field. 

An important outcome from [3[ is the observation that 
the problem of finding self-consistent IPs with single 
modulation can be reduced to finding the ground state 
solution in a (1+1)-D model. The relevance of this ob- 
servation is the understanding that the pairing dynam- 
ics of single-modulated solutions is essentially (1+1)-D. 
Thanks to this effective dimensional reduction, one can 
use known results from (1+1)-D theories like the Gross- 
Neveu (GN) QUI] or the NJL 2 HH, to tackle the prob- 



lem of finding self-consistent phases with 1-D modulation 
in a (3+l)-D NJL model. 

One might naively think that since the chiral spiral so- 
lution of the (1+1)-D QCD theory |l[ is directly related 
to the confinement of the interactions at the Fermi sur- 
face, one should not expect the same ground state in a 
(3+l)-D NJL model where confinement is absent. How- 
ever, as will become clear later, this argument is not valid 
because the dimensional reduction of the pairing dynam- 
ics always connects the (3+l)-D NJL models to lower 
dimensional theories with asymptotic freedom. 

As we shall show below, once we extend the NJL model 
of Rcf. [3J to allow a tensor-tensor four-fermion interac- 
tion term that preserves the chiral symmetry of the the- 
ory, the quarkyonic chiral spiral solution becomes more 
energetically favored than the real, crystalline solution of 
[3j. At the end of the paper, we shall discuss, using sym- 
metry arguments, what is the underlying physics making 
the connection between our modified NJL model and the 
confined quarkyonic matter. 

A modified (3+l)-D NJL Model. We start with 
the following two-flavor chiral invariant NJL Lagrangian 
in (3+1) dimensions 



c = i> (y* c*^ + Mmo) </> + g [(^/0 2 
+c [(v^v) 2 + (^r 75f 7^) 2 ] 



(ipiTj 5 ip) 2 



(1) 



where r are Pauli matrices associated with the SU(2) fla- 
vor symmetry, a^ v = |[7 a 'j7 1 ']) 7 m is taken in the chiral 
representation, and ip is the ANfN c dimensional quark 
spinor with Nr = 2. 



We consider the large N c limit, thereby avoiding the 
no-go theorem issue after the problem reduces to symme- 
try breaking in a lower dimensional theory [18|, ll9( . For 
simplicity, we are going to assume G = G' . We are in- 
terested in 1-D modulated condensates and without lost 
of generality can take the modulation in the z-direction. 
We shall consider the condensates {i>ip) = S(z) and 
(-ijja 03 tp) = D(z), because once the scalar condensate 
S(z) selects a preferential direction, nothing prevent in 
principle the D(z) condensate to exist. On the other 
hand, we can ignore the condensation of the second 
structures on each terms, because the pion condensate 
is known not to be favored [3j and a similar result can be 
shown for the ipir^a^ip structure. 

In the mean-field approximation the thermodynamical 
potential is 

™^^TrlnI[. 
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acts on Dirac and coordinate space. Using a proper regu- 
larization, one can express the thermodynamic potential 
in terms of the eigenvalues E n of Hmf as 
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The gap equation is obtained by minimizing the ther- 
modynamic potential with respect to M*(z). It can be 
expressed in terms of the normalized eigenvectors tp n (x) 
of Hmf for the eigenvalues E n [3(] , 
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(5) 
Following Rcf. [3j , we can take advantage of the trans- 
lational invariance in the (x,y) plane, which implies the 
conservation of the momentum component P± in that 
plane. The eigenfunctions of H M f can be labelled by 
the energy eigenvalues and the perpendicular momen- 
tum, V'A.p.L. • With the help of a Lorentz boost they can be 
expressed in terms of the eigenfunctions -0a, o for P_i_ = 0. 
This allows to construct the eigenvalue spectrum from 
the subspace spanned by {V'a.o} (see [3j for details). In 
this subspace the Hamiltonian ^ takes the form 
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which can be block-diagonalized to 
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where 



H 1D (M*(z)) 
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\ M{z) id z 



(7) 
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In the superconductivity literature Hmf is known as the 
Bogoliubov-de Gennes Hamiltonian [20( • It is the Hamil- 
tonian of the NJL2 theory whose solutions were exten- 
sively studied in [17] . Then, our problem has been re- 
duced to a direct product of two identical NJL2 Hamil- 
tonians, and hence we can identify the gap A in Ref. 
17 with M*(z) and take advantage of the results found 
there. The gap equation of the theory ([5]) can be solved 
using a resolvent approach [17], [2l| . Its general solution 
is given by a periodic array of chirally twisted kinks [l7j 



M*{z) = -\e 2lqz A 



a(\Az + iK'-id/2) 



a(\Az + iK')a(i9/2) 
x exp[i\Az(-i((i9/2) + ins(i9/2)) + i9 m /2] (9) 

where A = A(9,v) = -2isc(i6/4; v)nd(i6/4; v); sc = 
sn/cn, nd = 1/dn arc Jacobi elliptic functions, and a 
and £ are the Weierstrass sigma and zcta functions, cho- 
sen to have real and imaginary half periods: wi = K(V), 
and UI3 = iK'(v) = iK(l — v), with K(f) the complete 
elliptic integral. Both periods are controlled by a single 
(real) elliptic parameter < v < 1. The parameter A 
sets the overall scale of the condensate, and 1/A sets the 
length scale of the crystal. Hence, the general solution is 
specified by four parameters: a scale A, a phase param- 
eter q, an angular parameter 9 (9 e [0, 4K']), and the 
elliptic parameter v. 

For the extreme values of 9, the twisted kink solution 
reduces to a chiral spiral 



M = \e~ 2lqz 



(10) 



At zero temperature the energetically favored solution is 
the spiral with q = /1 [17] . This is the same quarkyonic 
chiral spiral solution found in quarkyonic matter at zero 
T Q. This result indicates that the (3+l)-D NJL model 
([T]) encompasses the same physics as the one taking place 
in a patch of the Fermi surface in the quarkyonic matter. 
This is encouraging, as one may then use this NJL the- 
ory to explore other regions of the parameter space that 
would be rather difficult to describe with QCD effective 
models with confined gluon propagators. Here again one 
can take advantage of the connection with the NJL2 
model and the existing results for this lower dimensional 
theory [l7J |. 

Ginzburg-Landau Analysis. We are particularly 
interested in exploring the region near the chiral transi- 
tion curve to determine the form of the most energeti- 
cally favored solution there. With this aim, we can use 



the Ginzburg-Landau (GL) expansion of the thermody- 
namical potential in powers of the condensate and its 
derivatives, 

* G l(M*) - * GL (0) 

=2{a 2 |M*| 2 +a 3 Im(ArAf')+a 4 (|Ar| 4 
+ \M*'\ 2 ) + a 5 Im((M"* - 3\M*\ 2 M*)M') 
+ a 6 [2|M*| 6 + 8|Af* | 2 |Af'*| 2 + 2Re((Af'*) 2 Af 2 ) 



\M" 



+ ...} 



(11) 



where the coefficients a n (T, /i) are functions of T and /1 
given in [l7| . Their explicit expressions are not relevant 
for our analysis. Expression (fTTj) is just twice the NJL2 
GL functional found in [l7| because the (3+l)-D Hamil- 
tonian (|7|) is the direct product of two Hir>(M*(z)), so 
its GL functional is the sum of two NJL2 GL functional. 
At each given order of the GL expansion, one can min- 
imize the GL functional with respect to M(z) to find the 
GL equation for the condensate at that order. Expanding 
up to «3 we obtain 



M' 



z^AT 
a 3 



=> M* = Ae 



(12) 



which is again a chiral spiral solution. One can verify that 
the mimimun solution at each order of the expansion is 
always a chiral spiral [l7[. Therefore, the chiral spiral is 
energetically preferred over the real crystalline solution 
of [3j in the region close to the chiral transition. 

Notice that the GL functional (fTTj). being twice the 
functional of NJL2, corresponds to a lower dimensional 
theory that has asymptotic freedom and continuous chi- 
ral symmetry [l7( . These are the same symmetries of the 
(1+1)-D QCD theory that describes the confined inter- 
actions in the Fermi surface of the quarkyonic matter. 

Discussion. In this letter we have unveiled a connec- 
tion between quarkyonic matter and a modified (3+l)-D 
NJL model that has an attractive chiral-invariant tensor- 
tensor four-fcrmion interaction channel in addition to the 
conventional scalar and pseudo scalar terms. Condensa- 
tion in the extra channel is favored once the spatial sym- 
metry is spontaneously broken by the scalar condensate. 
For single-modulated condensates, the pairing dynamics 
of our (3+l)-D NJL model becomes (1+1)-D and is de- 
scribed by the NJL2 theory. In the case of quarkyonic 
matter, the relevant pairing dynamics is described by 
(1+1)-D QCD. Since the symmetries of these two lower 
dimensional theories are the same, it is logical to expect 
that they posses the same ground state. As a matter 
of fact, the ground state of our (3+l)-D NJL model is a 
crystalline chirally twisted kink, which becomes equal, at 
zero temperature, to the chiral spiral solution found in 
the patches of the Fermi surface in quarkyonic matter. 

It is now easy to understand why the most favored so- 
lution for the NJL model ([1]) differs from the one found 



in [3j . The reason is that these two N JL models are con- 
nected to different lower dimensional theories which have 
different symmetries. The (3+l)-D Hamiltonian of the 
present case is connected to the Chiral GN or NJL 2 the- 
ory. In contrast, the Hamiltonian of Ref. [3j is connected 
to the GN model. The chiral symmetry is continuous in 
the first, but it is discrete in the second. For the first 
model the Hamiltonian reduces to the direct product of 
two NJL2 Hamiltonians, so the GL expansion is twice 
the NJL2 GL expansion. In the second case, the Hamil- 
tonian reduces to the direct product of H\r>{M{z)) and 
H\d{M*{z)), so the terms with odd-labelled coefficients 
in the GL expansion cancel out. Without these terms 
the minimum solution is never a chiral spiral, but a real 
function as shown in [3J. 

The connection of the proposed NJL model with the 
quarkyonic matter is rooted on the dimensional reduc- 
tion of the dynamics of fermions in the Fermi surface. 
The fermions close to the Fermi surface are the relevant 
degrees of freedom to produce the particle-hole conden- 
sate that yields the inhomogencous phase. Thus, de- 
spite the apparent differences between the initial (3+1)- 
D quarkyonic matter theory and the proposed (3+l)-D 
NJL model, one with confinement and the other without 
it, each of them reduces to a lower dimensional theory 
with confinement and the same continuous chiral sym- 
metry, thereby producing the same ground state. 

The proposed NJL model can be useful to explore the 
region in the fi — T plane occupied by the inhomogeneous 
phase, the kind of transitions that limit this region, and 
the possible existence and location of the QCD critical 
point. 

Acknowledgments: The work of EJF and VI has 
been supported in part by DOE Nuclear Theory grant 
Ide-sc0002179l 



[1] T. Kojo, Y. Hidaka, L. McLerran, and R.D. Pisarski, 

Nucl. Phys A 843 (2010) 37. 
[2] L. McLerran and R.D. Pisarski, Nucl. Phys A 796 (2007) 



83. 
[3] D. Nickel, Phys. Rev. D 80 (2009) 074025; Phys. Rev. 

Lett. 103 (2009) 072301. 
[4] K. Fukushima and T. Hatsuda, Rept. Prog. 

Phys. 74 (2 011) 014001; K. Rajagopal and 

F. Wilczek, |arXiv:hep-ph/0011333] T. Schafer, 

|arXiv:hep-ph/0509068| M. BubaTIa7~Phys. Rep. 407 

(2005) 205. 
[5] Y. Aoki et al., Phys. Lett. B 643 (2006) 46; JHEP 0906 

(2009) 088; H. Abuki et. al, Phys. Rev. D, 78 (2008) 

034034. 
[6] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 

345; ibid. 124 (1961) 246. 
[7] D.V. Deryagin, D.Y. Grigoriev, and V. A. Rubakov, Int. 

J. Mod. Phys. A 7 (1992) 659; E. Shuster and D. T. Son, 

Nucl. Phys. B 573 (2000) 434; B.-Y. Park et. al, Phys. 

Rev. D 62 (2000) 034015. 
[8] R. Rapp, E. Shuryak, and I. Zahed, Phys. Rev. D 63 

(2001) 034008; E. Nakano and T. Tatsumi, Phys. Rev. D 

71 (2005) 114006; N.V. Gubina et. al, Phys. Rev. D 86 

(2012) 085011. 
[9] S.Carignano, D. Nickel, and M. Buballa, Phys. Rev. D 

82 (2010) 054009. 
[10] H. Abuki, D. Ishibashi, and K. Suzuki, Phys. Rev. D 85 

(2012) 074002. 
[11] S. Klevansky, Rev. Mod. Phys. 64 (1992) 649. 
[12] B. Feng, E.J. Ferrer and V. de la Incera, Nucl. Phys. 

B 853 (2011) 213; Phys. Lett. B 706 (201 1) 232; E.J. 

Ferrer and V. de la Incera. larXiv: 1208.51791 [nucl-th], 
[13] E.J. Ferrer, V. de la Incera, and A. Sanchez, Acta Phys. 

Polon. Supp. 5 (2012) 679. 
[14] O. Schnetz, M. Thies, and K. Urlichs, Ann. Phys. (N.Y.) 

314 (2004) 425; O. Schnetz, M. Thies, and K. Urlichs, 

Ann. Phys. (N.Y.) 321 (2006) 2604; 
[15] M. Thies, J. Phys. A 39 (2006) 12707. 
[16] C. Boehmer, M. Thies, and K. Urlichs, Phys. Rev. D 75 

(2007) 105017; G. Basar, G.V. Dunne, Phys. Rev. Lett. 

100 (2008) 200404; C. Boehmer and M. Thies, Phys. 

Rev. D 80 (2009) 125038; F. Correa, G.V. Dunne, and 

M. S. Plyushchay, Annals Phys. 324 (2009) 2522. 
[17] G. Basar, G.V. Dunne, and M. Thies, Phys. Rev. D 79 

(2009) 105012. 
[18] N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17 (1966) 

1133; S. Coleman, Commun. Math. Phys. 31 (1973) 259. 
[19] E. Witten, Nucl. Phys. B 145 (1978) 110. 
[20] P. G. de Gennes, Superconductivity of Metals and Alloys 

(Addison- Wesley, Redwood City, California, 1989). 
[21] J. Feinberg and A. Zee, Phys. Rev. D 56 (1997) 5050. 



